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We show how the density-matrix numerical renormalization group (DM-NRG) method can be 
used in combination with non-Abelian symmetries such as SU(JV), the decomposition of the direct 
product of two irreducible representations requires the use of a so-called outer multiplicity label. 
We apply this scheme to the SU(3) symmetrical Anderson model, for which we analyze the finite 
size spectrum, determine local fermionic, spin, superconducting, and trion spectral functions, and 
also compute the temperature dependence of the conductance. Our calculations reveal a rich Fermi 
liquid structure. 
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I. INTRODUCTION 

Quantum impurity models, describing a quantum sys- 
tem with a small number of discrete states, the impurity, 
coupled to a continuous bath of fermionic or bosonic ex- 
citations, arise in a variety of contexts. A particularly 
important example is the Anderson impurity model, 1 rel- 
evant for describing magnetic moments in metals, trans- 
port through quantum dots, and for the treatment of 
correlated lattice models using dynamical mean field the- 
ory. While the standard version of this model has SU(2) 
spin symmetry in the absence of a magnetic field, gen- 
eralizations to settings with higher symmetry have also 
been studied. The SU(A) generalization of the Ander- 
son model emerged first in the context of heavy Fermion 
systems, 2 where large N expansions proved to be an ef- 
ficient way to model and describe magnetic atoms with 
orbital degeneracy. Studying these models in detail is not 
only useful in the context of heavy fermion systems, but 
it also represents the first step to understand the behavior 
of correlated cold atomic gases with SU(iV) symmetrical 
interactions. 3,4 

The SU(A) Anderson model can also be realized in a 
controlled way. In particular, the SU(4) model has been 
realized in various mesoscopic structures including car- 
bon nanotubes, 5,6 vertical quantum dots, 7 and more re- 
cently in the originally proposed double dot structures. 8 
Similarly, the SU(3) Anderson model could also be real- 
ized with quantum dot structures, though the proposed 
set-up is maybe somewhat more complicated. 9 

The SU(iV) Anderson model is defined in terms of TV 
local orbitals embedded in a conduction electron sea. Its 



Hamiltonian can be written in terms of the corresponding 
creation operators, S a (a = 1, ...,N) and the number 
operator h = d} a d a as 



H = e d h+ ^n(h - 1) + V ^(4^ a (0) + h.c.) 



chan 
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Here £d and U denote the position of the local orbital 
and the strength of interaction on it, respectively, and 
the level eft hybridizes locally with the fermions at its po- 
sition, destroyed by ^> a (0). The last term of the Hamilto- 
nian describes the kinetic energy of the conduction elec- 
trons. It generates the dynamics of the field ip a (0), and 
amounts in a broadening of the "atomic" level, Sd- 

In the present paper we show on the prototypical ex- 
ample of the SU(3) symmetrical Anderson model how 
the numerical renormalization group 10,11 (NRG) of Ken 
Wilson, one of the most versatile and reliable tools for 
treating quantum impurity models, can be adapted to 
fully take advantage of non-Abelian symmetries to re- 
duce computational costs. Within Wilson's procedure, 
one rewrites (1) as Hamiltonian of a semi-infinite chain, 
and diagonalizes it iteratively. 10 Using symmetries in the 
course of this diagonalization procedure is crucial: it al- 
lows computer memory to be used efficiently and en- 
ables one to reach the required numerical accuracy on 
relatively standard computers with reasonable runtimes. 
Eq. (1) obviously possess a SU(JV) x U(l) symmetry cor- 
responding to rotations in spin space and overall charge 
conservation. Here, we shall focus on the N = 3 case, 
classify states and observables while exploiting these 
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symmetries, and determine the spectral functions of sev- 
eral local observables. 

In an earlier work, a general framework has been set up 
and implemented to handle an arbitrary number of non- 
Abelian symmetries dynamically. 12 This formulation al- 
lowed us to build an open-access flexible density-matrix 
NRG (DM-NRG) code. 13 However, in Ref. 12 we con- 
sidered only combinations of certain rather simple sym- 
metries such as charge and spin SU(2) symmetries, Z 2 
or U(l) symmetries. In a group theoretical sense, these 
are simpler than SU(7V > 2) and some other discrete or 
Lie groups. For SU(2), irreducible representations (ir- 
reps) are labeled by the size of the spin. When " adding" 
two SU(2) spins, say Si and S 2 , each possible total spin 
S3 that satisfies the angular momentum addition rule 
15*1 — S 2 1 < S3 < S1+S2 can be obtained in precisely one 
way. More technically, in the decomposition of the direct 
product of two SU(2) irreps labeled by spins Si and S2 
into a direct sum of irreps, the number of times ns 1 s 2 ;S 3 
that any irrep labeled by spin S3 occurs in the direct 
sum, the so-called "outer multiplicity", is either or 1. 
This is, however, not true for SU(A > 2): the decompo- 
sition of the direct product Ti ® T 2 of two SU(iV) irreps 
into a direct sum can contain irreps with outer multiplic- 
ity nY 1 r 2 ;r 3 larger than 1, in other words, there may be 
n rir 2 ;r 3 incquivalent ways to construct the irrep from 
Ti and T 2 . Correspondingly, the Clebsch-Gordan coeffi- 
cients have a more complicated structure in this case, and 
the Wigncr-Eckart theorem, extensively used in the DM- 
NRG calculations, becomes also somewhat more compli- 
cated. Here we show how a general framework can be 
constructed to deal with this case, 14 and demonstrate it 
on the specific example of the SU(3) Anderson model. 

We note that another general approach towards ex- 
ploiting non-Abelian symmetries such as SU(iV) or 
Sp(iV) within NRG, and more generally for tensor net- 
work methods, has recently been published by A. Wc- 
ichsclbaum. 15 It is formulated in the language of matrix 
product states, and may be regarded as complementary 
to our own, which is phrased within the more traditional 
formulation of NRG. We emphasize, though, that both 
approaches are fully equivalent, in that precisely the same 
NRG assumptions and approximations are made in both; 
they differ only in the data structures used for internal 
book-keeping in the numerical codes. Their relation is 
briefly scotched in App. D. 

Our DM-NRG calculations require explicit knowledge 
of Clebsch-Gordan coefficients. Whereas these are known 
in closed form for SU(2), this is not the case for SU(A > 
2). However, an efficient numerical algorithm for their 
evaluation has recently been developed, 16 which we use 
here. 

This paper is structured as follows: In Sec. II we out- 
line our approach for exploiting non-Abelian symmetries 
in DM-NRG calculations. In Sec. Ill we apply it to the 
SU(3) Anderson model; in particular, we present results 
for the conductance through a quantum dot described by 
this model, and for various local spectral functions. Our 



conclusions in Sec. IV arc followed by four appendices, 
that summarize some basic facts of SU( N) representation 
theory, and some recursion formulas involving Clebsch- 
Gordan coefficients, respectively. 

II. DM-NRG WITH NON-ABELIAN 
SYMMETRIES 

As stated in the introduction, the formalism presented 
in Ref. 12 applies only for a special (though relatively 
large) class of symmetries. Therefore, let us review here 
the most important formulas and the structure of the 
NRG calculations for the more general case, where so- 
called outer multiplicities are also considered. 1 



A. Local symmetries on the Wilson chain 

Let us start by first discussing the general structure 
of the symmetries of the Wilson chain. The first step in 
Wilson's procedure of solving a quantum impurity prob- 
lem is to perform a Gram-Schmidt orthogonalization and 
rewrite the Hamiltonian in a "tridiagonal" form, 11 

00 

+ Hn+i). (2) 

Here Hq contains the local, interacting part of the Hamil- 
tonian, while the rest of the chain represents the conduc- 
tion electron (bath) degrees of freedom, coupled to it. 
The on-site terms H. n +i are many times missing; they 
typically appear for more sophisticated electronic densi- 
ties of states and can also account for superconducting 
correlations. In our case, one could take, e.g., 

Ho =s d h+ ^n(n - 1) + 1/]T(4/I 01 + h.c.) , (3) 

a 

as Hq, with f$ ~ ip a (0) a properly normalized on-site 
fcrmion, and the hopping terms r ntn +i would read, 

r„,„ +1 =^/ l W( / Ht / [n+i] +h . c .) (4) 

a 

Here /<L n+1 ' annihilates a fermion of SU(3) spin a at site 
[n], and the hopping amplitudes h<- n ' decay exponentially 
along chain, thereby leading to energy scale separation. 
Eq. (2) is then diagonalized iteratively using the recursive 
relation 

H n = H n _i + T„_i >n + H n . (5) 

In the following, we shall assume that H (and H n ) 
are invariant under the direct product of ns symmetry 
groups, 

G = Qi x g 2 x ■ ■ • x g ns . (6) 
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This means that for any group element g\ £ S\, with A = 
1, . . . , ng, there exists a unitary operator U\(g\) on the 
Fock space, leaving H invariant. Here we do not need to 
make much restriction on the groups: our considerations 
hold for any group which acts on the chain locally at each 
lattice site n 

n s n s 

u{g) = n ^a(.9a) = n nx«G?A) , (7) 

A=l A=l n 

and for which the Wigner-Eckart theorem holds. 18 Given 
the above group structure, we can group all eigenstates 
and also operators into multiplets, each of which trans- 
forms according to a certain representation of Q, 

r = {r 1 , r 2 , r ns } o r 1 ® r 2 ® ■ ■ ■ ® r ns . (8) 

The "quantum numbers" T A , which label the various ir- 
reducible representations (irreps) occurring in T, can be 
spin labels, charges, or can label different irreps of some 
point group. States within a multiplet are then labeled 
by a set of internal indices, 7 = {7 1 , 7 2 , j ns }, with 
the internal labels running from 1 < j x < dim(r A ). A 
given multiplet i that transforms according to the repre- 
sentation Ti thus consists of 

n s 

dim(i) = dim(r. i ) = J| dim(r\ A ) (9) 

A=l 

degenerate states. 

States belonging to a product of two representations, 
r A © , often appear in the calculations. As outlined in 
the introduction, similar to spins, these can be decom- 
posed into irreps. However, one irrep may appear several 
times in this decomposition, 

r A ©r A -> ...©r A ©...©r A ©... . (10) 

times 

Accordingly, in the most general case, the resulting states 
r A <£> Tp — >• r A must be labeled by a so-called " outer mul- 
tiplicity" label, a a = l,..,n A p . q . Correspondingly, prop- 
erly transforming multiplets may be constructed from the 
product of two multiplets as 

|r,7>a= ^2 (ri,7i; r 2,72|r,7) a |ri, 71) ® |r 2 , 72) , 

Ai.A 2 

where a = {a\} denotes the composite multiplicity la- 
bel, and the generalized Clebsch-Gordan coefhcients are 
defined as 

ns 

(ri, 7 i;r 2 ,72 1 r, 7 ) Q = J] (r A ,7i A ;r 2 A ,7 2 I r A , 7 A ) QA ■ 

A=l 

The outer multiplicity also appears in the Wigner- 
Eckart theorem. The latter states that the matrix ele- 
ments of an operator multiplet, i.e. a set of operators 



{Ar A , 1A } transforming under transformations U(g) as a 
multiplet r^, are determined almost entirely by repre- 
sentation theory, and can be expressed in terms of the 
Clebsch-Gordan coefficients as 

(i,T u7i \A TAnA \j,T j , 7j ) = 

^2(r A , 7A ;r j>7j \ r i>7i y a (i\\ a\\ j) a . (13) 

OL 

Here multiplets i and j transform according to the irreps, 
Fj and F,. Thus, according to the Wigner-Eckart theo- 
rem, all matrix elements are determined by only a few 
reduced matrix elements, (i \\ A \\ j) a , labeled just by 
the outer multiplicity labels, a characterizing how many 
times the representation appears in the product of Tj 
and r^. For many commonly used symmetries as SU(2), 
e.g., the outer multiplicity is just always one and the la- 
bel a can be dropped. However, it is needed for, e.g., 
SU(N > 3) or even for cubic point groups. 

B. Wilson's NRG with symmetries 

In course of the NRG procedure, one diagonalizes 
Eq. (2) iteratively. The eigenstates of the Hamiltonian 
H n of a chain of length n can be grouped into multiplets, 
with each multiplet i transforming according to a cer- 
tain representation Ti = {T 1 , .., r" s }. Having computed 
the approximate eigenstates (block states) \i, Ti, 7.;, 
of H n _i, one proceeds to construct eigenstates of H n . 
To do that, one first appends to the chain the multi- 
plets {\n, T l ° c , 7^ oc )}, spanning local Hilbcrt space at site 
n, and then constructs properly transforming multiplets 
{|it, r„, 'ju^al } by making use of the Clebsch-Gordan co- 
efficients, Eq. (11): 

| A i,r^, 7 ^ oc ) ® |i,r i ,7i> [n - 1] -> \u,T u , lu )W . (14) 

Notice that a new multiplet u now also carries an outer 
multiplicity label, a u : This specifies the representation 
according to which T u has been produced from and T a . 
The advantage of using these states is that H n is diagonal 
both in r„ and in the internal labels, j u . Therefore, it is 
sufficient to compute only the corresponding irreducible 

[n] 

matrix elements (u \\ H n || v) in each symmetry sector 
(block) separately, and diagonalize H n sectorwise by a 
unitary transformation to obtain the corresponding new 
eigenstates, 

|u J r , 7u )W^|? ) r s , 7 5>w. (15) 

As explained in App. D, this iterative procedure leads to 
a matrix product state (MPS) with a peculiar structure, 
reflecting the symmetry of the Hamiltonian. 

The most difficult part in the procedure above is to 

[nl 

determine the matrix elements (u || H n j| v) . These 
can be constructed by noticing that each state u and v 
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has been constructed from the eigcnstates of H n _i and 
% n , i,fi—tu and j, v — > v, and therefore 

(u || H n || «)W = S^iE?- 1 + «£) + (u || r„_ 1>n || «)["!, 

with e" being the eigenenergy of ~H n . The matrix ele- 
ments of T„_i jn can be worked out by assuming that the 
hopping part consists of some fermionic or bosonic cre- 
ation operators C^\. , transforming again according 
to some irreps T a , 



T~n—l,n 



7, 



(16) 

Here a labels the different "hopping operators", and 
ha ^ the corresponding hopping amplitudes. Notice 
the somewhat unusual way this hopping term is written: 
C-aXa-ia / is a "creation operator", which transforms 
according to the representation T a , while (C r „)' O / 



is an " annihilation operator" , transforming according to 
the conjugate representation, T a *. We remark that for 
charge SU(2) symmetry e.g., the "creation" operator 
multiplet C Qi r a ,7 a is a Nambu spinor, and contains both 
/ and f> operators. 20 The number of hopping operators 
may depend on the symmetry used: for a chain of spin 
1/2 fermions treated in terms of SU(2) x U(l) symmetry, 
e.g., one has a single hopping operator of spin 1/2 and 
charge 1, while if only the charge symmetry is used then 
one has two hopping operators of charge 1 , corresponding 
to the spin up and spin down directions. In our example 
of the SU(3) Anderson model we have a single hopping 

operator, and C^ «-> {/<L™^}- Assuming then that the 
reduced matrix elements of the creation operators act- 
ing on site n — 1 of the chain, (u || C a n 1 ' || v)a , 
and those of the local creation operators at the added 
site, (v || Ca II M}/3 are known, one can use the Wigner- 
Eckart theorem to express (u \\ r n -i,n II as 



(u || T n -i,n || v) n = S r , 



(i II || j)^ 11 (v II || tf p D (a, a, ft u, v) + h.c. 



(u «- Df- . (17) 



Here the outer multiplicity labels a (,5) label inequiva- 
lent ways in which I\ (r„) appear in the product r a <S>Tj 
(r a (g) r^). The coefficients D (a,a, (3;u,v) can be ex- 
pressed in terms of Clcbsch-Gordan coefficients, and are 
given by Eq. (C3) in Appendix C. Similar expressions 
hold for the matrix elements of "block operators" (A), 
acting somewhere on the first n — 1 sites of the chain, 
and those of "local operators" A^ n \ acting only on the 
last site of the chain. 19 For a block operator we have, 
e.g., 

(« II A || v)f = £(i || A || F(a,P;u,v) <W , 

a 

(18) 

while for the local operators the following equation holds, 

(u || AW |j w>W = ]T> II A [n] || ^} tt K(a,/3;u,v) S itj , 

a. 

(19) 

with the coefficients -F (a, /?; u, v) and -fC (a, j3; u, v) given 
in Appendix C. Here again, the outer multiplicity la- 
bels j3 label inequivalent ways in which T u appears in 
Ta ® T v , while a labels similarly inequivalent ways how 
I\ (Tfj,) can be constructed from Ta and Fj (Fa and r„). 
Similar to D, the coefficients F and if are again deter- 
mined only by symmetry, and can be expressed in terms 
of Clebsch-Gordan coefficients. As a last step of the it- 
eration, the reduced matrix elements (u || A^ || v)^ 
and (u || A \\ v)$ need be transformed to the new basis, 



|z,r^,7^}. This is performed by using precisely the same 
unitary block transformations as the ones used to diag- 
onalize the Hamiltonian H ni (15), without affecting the 
outer multiplicity labels /3. 

Wilson's diagonalization procedure can be carried out 
then based upon the equations above: In a given itera- 
tion, one takes the lowest lying states of iteration [n — 1] 

and their matrix elements (i || C a n ^ || j)g , computes 

from these and from the matrix elements (/x || cl™' \\ v)p 
the Hamiltonian (u \\ H n || v)^. Then diagonalizing 
(u || H n || one obtains low-lying eigenstates of H n 

and determines their matrix elements (i | cl™' || 
Continuing this procedure for larger and larger values 
n, one obtains better and better approximations for the 
ground state of H = H n ->oo and the lowest lying eigen- 
states. 



C. FDM-NRG approach 

So far, we discussed essentially Wilson's original NRG 
approach in case of general symmetries. In practice, how- 
ever, one often needs to go beyond Wilson's RG and 
use the so called DM-NRG method, 21 whereby a forward 
NRG run is first performed to obtain the density matrix 
(DM) of the system, and then a backward NRG run is 
made to compute physical obscrvablcs. Moreover, to sat- 
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isfy spectral sum rules, a complete basis set 25,26 has to 
be used, as first implemented in the context of DM- NRG 
in Refs. 22 and 23. In the full density-matrix NRG ap- 
proach (FDM-NRG) of Rcf. 23, the full density matrix of 
the entire chain is expressed in the complete basis, which 
yields an improved treatment of finitc-tcmpcraturc prop- 
erties. Let us now briefly discuss how symmetries can be 
implemented in the FDM-NRG approach. (For a compli- 
mentary formulation of the same strategy using matrix 
product states, see Ref. 15 .) 

First, to satisfy the necessary completeness relations, 
we consider a chain of N sites and introduce "environ- 
ment" states e for a each state discarded in iteration [n] 
(i G £), 22 > 23 



(20) 



Here the states e„ form an orthonormal basis for the 
remaining N — n sites of the chain, and their internal 
structure is irrelevant for the remaining discussion. 

The states (20) form a complete basis on the Wilson 
chain, 25 ' 26 and can be used to construct the density op- 
erator as follows 23 



N 

e = E * [B] , 

71=0 

* [Bl - E E 



(21) 



"1^5 Ti, 7i ) ^ ^ ^ (j'j , 7z , C n | . 



set of truncated reduced density matrices, 



R [n] ee Tr 



(23) 



By symmetry, the reduced density matrices are invariant 
under the symmetries used, and have therefore a block- 
diagonal structure in the representation indices. 12 The 
matrices 

R[n-1] 

can be constructed iteratively from W n > 
by tracing out site n and then adding the contribution of 
states discarded in iteration (n — 1) — > n. The contribu- 
tion of the kept states (K) reads: 



(i || || j) 



dim(u) 
dim(i) 



(u || || w)W . (24) 



Here the tilde indicates that the summation runs over 
states u and v having the same symmetry, and con- 
structed from states i and j by adding the same local 
state, i ® — y u and j ® ji ^ v. The subscript indicates 
that i and j are both kept states . The discarded piece 
of is then simply 



(25) 



To gain insight to the dynamical properties of a quan- 
tum impurity, one usually computes the retarded Green's 
functions for some operator multiplcts A 1A and B lB , 



Here /3 = 1/ksT and the partition function is expressed 
as 

N 

Z = ^^dim(i)e-^"^-™ ) (22) 

n=0 ieD 

with d the dimension of the Hilbert space at each added 
site of the Wilson chain, and d N ~ n the dimension of the 
space of the "environment" states, e n . We remark that in 
the last iteration all states are considered to be discarded, 
while in the first few iterations there are typically no 
discarded states yet. 

To compute local observables and spectral functions of 
observables at the impurity site, one traces out step by 
step the environment states, and introduces the following 



G£, 7 » = -- J Tr{ [(A 7 jt( t))j B 7B (o)] J e^dt , 

(26) 

with £ = — (£ = +) referring to commutators (anti- 
commutators) appearing for bosonic (fermionic) opera- 
tors. By symmetry, to have a non-vanishing value, A 
and B must both transform according the same repre- 
sentation, = Tb, and 7 a = 7b must also be satisfied, 
G^ t 7B = 5 1AnB G r j$ B . Notice that Eq. (26) is defined in 

terms of (A 1A )\ transforming according to the conjugate 
representation, rj\ 24 

The expression above can be evaluated in terms of the 
truncated density matrices, i?^ , and the reduced matrix 

elements (i || A |j j}a and (i || B || j)a of the operators 
A and B, to obtain 12 ' 26 
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< -Mt,s( 2 ) 



N 

E 

n=0 i£D,K (J,k) 



E 



E ( t 'll ij[nl Hi> [n] x 

(K,K) 



E 



(fc || A ||i)L n] *(fc || B || i)L™ ] dim(fc) 



El dim(A) 



(j || B || fcg(j M||fc)k" J * dim(i) 
z - \{E? + E™) + El dim( A) 



(27) 



This expression provides an efficient way to compute 
spectral functions. Notice that it contains only the re- 
duced matrix elements and the dimensions of the multi- 
plets involved. 



III. STUDY OF THE SU(3) ANDERSON MODEL 

To demonstrate how the scheme presented above 
works, we apply it to study the repulsive SU(3)- 
symmetrical Anderson model, defined already in the In- 
troduction. We perform our calculations for a conduc- 
tion band with a uniform local density of states be- 
tween energies W > e > —W with the bandwidth set 
to W = 1, and use the corresponding hopping ampli- 
tudes ftW ~ (1/2) (1 + A-^A"™/ 2 in (4). In this case, 
the width of the (nonintcracting) level is approximately 
given by A = irg c V 2 with g c = 1/2W the local density 
of states at site of the Wilson chain. 

As mentioned before, the Hamiltonians (1), (3) and (4) 
have a U(l) x SU(3) symmetry in the charge and flavor 
sectors, respectively. Correspondingly, multiplets of the 
Hamiltonian are characterized by a charge and a flavor 
quantum number. The charge quantum numbers Qi are 
simply identical to the total charge, 28 



Q 




dl. d n 



(28) 



conserved by (1), (3), and (4). 

Labeling the SU(3) representations and the states 
within an SU(3) multiplct is somewhat more compli- 
cated. The flavor quantum numbers Fi can be (and are 
usually) represented by Young tableaux, characterized 
by two non-negative integers in case of SU(3) (see Ap- 
pendix A). Young tableaux provide a nice pictorial way 
to multiply and decompose representations, or calculate 
their dimensions. For our numerical calculations, how- 
ever, we had to construct explicitly the basis states of 
SU(3) representations and to compute the corresponding 
Clebsch-Gordan coefficients. 29 This we carried out using 
the so-called Gelfand-Tsetlin patterns, briefly discussed 
in Appendix A. The states obtained this way are ana- 
logues of the canonical SU(2) basis states, created by 
the raising and lowering spin operators, S^. Gelfand- 
Tsetlin patterns are in one-to-one correspondence with 
the Young tableaux, but they allow for a simpler explicit 
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FIG. 1. (Color online) The occupation number (upper panel) 
and the conductance (lower panel) as function of e«j and for 
different broadening parameters A. The temperature is T = 
and the Coulomb interaction is fixed to U — 0.2 W. 



construction of the basis states. For more details, we 
refer the reader to Ref. 16. 

Similar to the cigenstates and multiplets, we also need 
to group operators into SU(3) multiplets and character- 
ize them by appropriate SU(3) quantum numbers. This 
classification of the most important operators is sum- 
marized in Table I. In terms of SU(3), there is only a 

single hopping operator, cJy™' O /Jp^, which transforms 
according to the defining SU(3) representation and has 
charge Q = 1, similar to the creation operator of the lo- 
calized level, S a . From d) a we can also construct various 
local operators of interest. The spin operators, ~ d) a dp 
form, e.g., an 8-dimensional charge Q = operator mul- 
tiplet in terms of U(l) x SU(3), while the charge Q = 2 
pairing operators ~ dj^dp transform according to a 3- 
dimensional SU(3) representation. Finally, the "trion" 
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Operator 



JO! ) U Q 



"a u /3 



(Q,F) 



(!.□) 



444 ( 3 >") 1x1 



Dim 



1 x 3 



1 x 3 



1 x 8 



Components 



CD -> ft 



-»■ 44 



44 



44 



1 








1 




J2_ 




1 




j5_ 





4^3 



-4^2 



1 




_3_ 





4di 



^ ^ 75 - 4d2 



+ 



24d 3 ) 



5-1 -> -4^2 

3 4rfi 



• — > 444 



TABLE I. Irreducible tensor operators for the SU(3) symme- 
try. 



operator 444 nas cnar g e Q = 3 and is an SU(3) sin- 
glet. 



A. Numerical results 

In our runs we have kept about 250 multiplets at each 
iteration, corresponding to approximately 1500 states on 
average, while the Wilson parameter was fixed to A = 2. 

In Fig. 1(a) we present the occupation of the localized 
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FIG. 2. NRG finite size spectrum at half-filling and T=0. 
Upper (lower) panel represents the even (odd) part of the 
spectrum. In both plots we represent the lowest 50 energy 
levels. The parameters are: A = 2, U = 0.2 If, £d = — U and 
U/ivA — 5. The convergence was reached after approximately 
27 iterations. 



level as a function of the energy Ed- The localized level 
can accommodate up to 3 fcrmions. For large enough 
U/A's, these fermions enter the local level one by one as 
£d is decreased, and a Coulomb staircase is clearly visible. 
In the range — U < £d < there is approximately one 
electron on the level. The SU(3) spin of this electron 
is screened by the conduction electrons, and an SU(3) 
Kondo state is formed. A similar, hole-like SU(3)* state 
emerges for — 2U < Ed < —U. 

The most intriguing region is though the mixed va- 
lence region, Ed ~ —U. For Ed — —U the ground state 
of the isolated impurity (r = 0) would be sixfold degen- 
erate due to electron-hole symmetry, connecting the two 
SU(3) triplets, {<4|0)} and {d+,d^|0}}. Valence fluctua- 
tions produce a state with all these states strongly mixed 
by quantum fluctuations. Fig. 2 and Table II show the 
flow diagram of the NRG levels, and the SU(3) classifica- 
tion of the asymptotic finite size spectrum, respectively. 
A detailed analysis reveals that this finite size spectrum 
can simply be understood as the finite size spectrum of 
three chiral fermions with a phase shift S = tt/2 at the 
Fermi energy. This phase shift is indeed in full agreement 
with the Friedel sum rule, 35/tt = (n), and the occupa- 
tion (n) = 3/2 assured by electron- hole symmetry. 

Similarly, the SU(3) Kondo spectrum, displayed in Ta- 
ble III can be understood as the finite size spectrum of 
three chiral fermions with a phase shift S m 7t/3, implied 
by the Friedel sum rule and the occupation (n) w 1. 

The crossover between the two SU(3) Kondo regimes 
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TABLE II. Energy spectrum for the electron-hole symmetric, 
mixed valence point. El stands for the finite size energy scale, 
E L = 2-k/L = 1.311. Qo = 0(3/2) for even (odd) iterations. 



through the mixed valence regime is maybe best captured 
by the local level's spectral function, shown in Fig. 3. In 
the SU(3) and SU(3)* Kondo regimes we find a narrow 
Kondo resonance of an exponentially small width pinned 
somewhat asymmetrically to the Fermi energy, u> = (see 
inset), and two Hubbard peaks. At the mixed valence 
point, £d = —U on the other hand, a relatively broad 
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TABLE III. Energy spectrum in the Kondo regime. Qo 
(— 3/2) for even (odd) iterations. 



and symmetrical resonance of width ~ T appears at the 
Fermi energy, and the charging peaks are absent. 



B. Conductance 

Let us now consider the mesoscopic structure proposed 
in Rcf. 9, sketched in Fig. 4. As explained in Ref. 9, this 
quantum dot structure would possibly realize a Kondo 
state with an approximate SU(3)-symmetry. The con- 
ductance through the device can be directly related to 
the spectral functions of the e?-level, and for a symmetri- 
cal device one finds 



(29) 



Here A(lj) is the spectral function of the d a operator 



of the localized level, A(uj) = 



'd a d 



(jjj)/tt, with 



G d d t (w) the Green's function defined in Eq. (26). The 

corresponding T = temperature linear conductance is 
shown in Fig. 1(b) as a function of the energy e^. The 
conductance reaches its maximal value in the mixed va- 
lence regime, £d ~ — U, and displays SU(3) Kondo effect 
related plateaus in the regions ~2U < Ed < —U and 




co/U 



FIG. 3. (Color online) a) The normalized spectral function 
for the on-site creation operator d} a for two different fillings 
n = 1.5 (half filling, red dashed) and n — 1 (1/3 filling, black 
solid). In the singly occupied case, n = 1, the dot is in 
the Kondo regime, and the resonance is shifted away from 
the Fermi energy, w = 0. The value of the spectral func- 
tion at the Fermi energy is determined by the Friedel sum 
rule, nAA(0) — s- 3/4. b) Evolution of the normalized spec- 
tral function as function of U/irA in the Kondo regime, for 
Ed/U — —0.5. The deviation from the Friedel sum rule is less 
than 1%. 



— U < e c i < 0. The T — temperature conductance val- 
ues observed can be understood in terms of the Friedel 
sum rule, relating the total occupation of the ef-level to 
the phase shift of the conduction electrons, which yields 
36/ tt = (n) for the SU(3) Anderson model. 9 In a Fermi 
liquid state, — at T = temperature, — the conduc- 
tance can be computed using the Landauer-Buttiker for- 
mula, and in the geometry of Fig. 4 is simply related 
to the phase shift as G(T = 0) = {e 2 /h) sin 2 ((5). This 
explains the value Gsu(3) ~ (3/4) e 2 /h observed in the 
Kondo states; there the occupancies are (n) « 1 and 
(n) « 2, corresponding to phase shifts 5 = ±7r/3, and 
the previously mentioned value of the conductance. At 
the mixed valence point, Ed = —U, on the other hand, 
we have (n) = 3/2, implying a phase shift 6 = tt/2, and 
a maximal conductance, Ggu(3) = e 2 /h . 

Increasing the temperature, the conductance is quickly 
suppressed in the Kondo regimes, and three Coulomb 
blockade conductance peaks emerge at the points of 
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FIG. 4. (Color online) Possible realization of SU(3) Kondo 
states using four edge states and three quantum dots. The up- 
per edge state is splitted and allows transport measurements 
between the two upper external leads. 9 

charge degeneracy, as shown in Fig. 5. The central peak 
corresponds to the transition n = 1 O 2 while the two 
side peaks correspond to charge fluctuations n = 1 -o- 
and n = 2 -H- 3, respectively. In Fig. 5 we also displaye4 
the temperature dependence of the conductance at the 
mixed valence point and in the SU(3) Kondo regime. The 
Kondo temperature of the SU(3) Kondo state is clearly 
much smaller than the mixed valence energy scale even 
for these moderate interactions. This Kondo tempera- 
ture can be estimated by first doing perturbation theory 
in V and performing a Schricffcr- Wolff transformation 
and then carrying out a renormalization group analysis. 
This analysis yields a Kondo temperature of 

T£ U(3) » D e-Va* , (30) 




FIG. 5. (Color online) Finite temperature conductance as 
function of Sd- Inset: The temperature dependence of the 
conductance on a logarithmic scale, for two different fillings, 
n = 1.5 (half filling, red dashed) and n = 1 (1/3 filling, black 
solid) . In both panels the Coulomb energy and the broadening 
to the leads were fixed to U = 0.2 W and U /tvA = 5. 
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FIG. 6. (Color online) Upper panel: Spectral function of the 
spin operator d^dp on a logarithmic scale for two different 
fillings, (n) = 1.5 (half filling, red dashed) and (n) ~ 1 (1/3 
filling, black solid) for U = 0.2 W . In the singly occupied 
case, (n) ~ 1, the dot is in the Kondo regime. The inset 
indicates the linear decay in As(lo) in the small frequency 
limit. The blue, dashed-dotted line is a guideline for the eye. 
Lower panel: The spectral function of the spin operator in 
the Kondo regime, (n) ~ 1, and for three different ratios 
U/wA: U/ttA = 15 (solid black line), U/ttA = 10 (dashed 
red line) and U/ttA = 5 (dashed-dotted green line). The inset 
indicates the universal scaling collapse of the spin spectral 
function in the Kondo regime. 



with the dimensionless coupling A expressed as 

A A 



A 



nE 4 



(31) 



in terms of the level width A and the "ionization en- 
ergies" E + = U + Ed and E- = — e^, and Dq ~ 
min(i? + , a high energy cut-off. 



C. Spectral functions of local operators 

The behavior of strongly correlated cold atomic and 
heavy fermion lattice systems can often be understood in 
terms of a self-consistent quantum impurity model (dy- 
namical mean field theory) . Within this picture, the local 
response functions of the quantum impurities may drive 



superconducting or magnetic phase transitions, or lead to 
even more exotic quantum phases. In this subsection, let 
us therefore analyze the spectral properties of the SU(3) 
Anderson model, and investigate the local spectral and 
response functions of its spin, pairing, and trion opera- 
tors. 

As discussed before, the spin operators, <v a dp trans- 
form according to an 8-dimensional SU(3) representation. 
As shown in Fig. 6, their spectral function displays Fermi 
liquid properties (see inset), and behaves very similarly to 
the spin spectral function of a standard SU(2) Anderson 
model. 30 In the mixed valence regime, for A <C U charge 
fluctuations to the state n — and n = 3 arc frozen out, 
and at low energies the only relevant energy scale is A; 
correspondingly, the spectral function exhibits a broad 
resonance at lj ~ A (extending up to u> <~ U), and decays 
linearly to zero for small frequencies, Ag~ 3 ^ 2 (cj) ~ lu/ A 2 . 
By Hilbert transform, this amounts in a spin susceptibil- 
ity: XS ~ 1/A. In the Kondo regime, (n) rj 1, on the 
other hand, two separate scales can be distinguished. Be- 
low lu ~ min([7, |e<j|) charge fluctuations are frozen and 
a clear Kondo anomaly can be observed as a logarithmic 
increase of the spectral function, as indicated by the ar- 
row in Fig. 6. The Fermi liquid behavior only emerges 
below the Kondo scale, u> < Tk "C A, U. In this Kondo 
regime A^ ond °(u) ~ oj/T^, and correspondingly, a sus- 
ceptibility xs ~ l/?if is found. 

The spin spectral function becomes universal in the 
Kondo limit, Tk <C A, U in the sense that the T = 
temperature dynamical susceptibility scales as xs^) = 
(1/Tjf) f(u/T K ,{n)), with f(u/T K ,(n)) a function, 
which only slightly depends on the occupation of the 
level, (n) 1. This is demonstrated in the lower panel 
of Fig. 6 for the imaginary part of the susceptibility, 
Xg(o>) = —71-^45(0;), computed for different values of 
U/irA. Numerically we define Tk as the half width at 
half maximum of A(uj), the spectral function of the d 
operator of the localized level. 
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0.05 




FIG. 7. (Color online) Spectral function of the supercon- 
ducting operator for U = 0.2 D and (n) = 1.5 (half filling, red 
dashed) and (n) « 1 (1/3 filling, black solid). 
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The correlations of the pairing operators, d^dp, be- 
have somewhat similar to those of the spin in the sense 
that at small frequencies a linear frequency dependence is 
found, corresponding to a Fermi liquid state with a con- 
stant pairings susceptibility (see inset of Fig. 7). How- 
ever, the local pairing operators lead out of the local low 
energy charging states both in the mixed valence and in 
the Kondo regimes. Therefore, the amplitudes of their 
spectral functions as well as the corresponding pairing 
susceptibilities obtained by a Hilbert transform, are typ- 
ically small, and they are not expected to drive any tran- 
sition. 



As a final example, we display the spectral function of 
the SU(3) singlet trion operator, = d\d\d\ m Fig- 8. 
The trion operator plays an important role in the attrac- 
tive case, 31 however, in this repulsive model it is a highly 
suppressed operator. In the Kondo regimes, it has a non- 
zero spectral function only because the Kondo states have 
a small (~ A/U) admixture of the empty and the triply 
occupied states, respectively. This explains why the am- 
plitude of the signal is in this case smaller in the mixed 
valence regime. It also explains the strong electron-hole 
asymmetry in the Kondo regimes. In the (n) ~ 1 regime, 
shown in Fig. 8, e.g., the admixture of the n = state 
is relatively large, ~ A/U, while the n = 3 state has 
a much smaller weight, ~ (A/U) 2 . As a consequence, 
most of the spectral weight appears on the particle-like 
side of the spectral function, u> > 0. At small frequencies 
the spectral function decays as ~ uj 2 . This is in agree- 
ment with Fermi liquid theory, which would predict a 
(T(t)TT(O)) <~ 1/t 3 decay of the trionic correlation func- 
tions at very long times. 
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FIG. 8. The spectral function for the trion operator for two 
different fillings n = 1.5 (half filling, red dashed) and n = 1 
(1/3 filling, black solid). The inset represents the same data 
on a logarithmic scale. The trionic spectral function scales as 
uj 2 in the small energy limit. 



IV. CONCLUSIONS 

In this paper, we showed how to extend the DM-NRG 
scheme of Rcf. 12 to symmetries with outer multiplici- 
ties. As an application, we performed a detailed DM- 
NRG study of the SU(3) symmetrical Anderson model by 
first incorporating SU(7V) symmetries 16 in the Open Ac- 
cess Budapest DM-NRG code, 13 and then performing the 
numerical calculations using the complete U(l) x SU(3) 
symmetry of the model. A similar extension has been 
carried out within the matrix product state (MPS) ap- 
proach parallel to this work. 15 

The properties of the SU(3) Anderson model do not 
differ so much from those of the original Anderson model. 
As also discussed in Refs. 9, for U > A four distinct 
charging regions appear: the featureless empty and fully 
occupied regions ((n) w 3 and (n) 0), and two Kondo 
regions of occupancies (n) sa 1 and (n) ~ 2, respec- 
tively. The two SU(3) Kondo regions behave similarly: 
they are characterized by phase shifts 8 ~ ±7r/3, as veri- 
fied from the finite size spectrum, and correspondingly, a 
Kondo resonance shifted away from the Fermi energy. In 
these Kondo regimes, the susceptibility has a universal 
form, x( w ) = fiu/Tx, (n))/Tf S -, with a scaling function 
f(x, (n)) very similar to the one emerging in the SU(2) 
Anderson and Kondo models. For completeness, we also 
studied the spectral properties of other local operators 
such as pairing or the trion operators. Both of them 
turn out to have a small amplitude for A <C U, and 
exhibit simple Fermi liquid properties below the mixed 
valence and Kondo scales, respectively. Therefore, away 
from half filling, a magnetic instability is expected to pre- 
vail on a lattice in the SU(3) Hubbard model, in general 
agreement with the results of Gutzwiller calculations at 
low temperatures. 2 

The SU(3) Kondo regions arc separated by a mixed 
valence state, which again has a Fermi liquid character 
with a Fermi liquid scale of the order of the level width, T. 
Here we find a phase shift 6 = 7r/2, in agreement with the 
expectations based upon the Friedel sum rule, but apart 
from that, and the emerging electron-hole symmetry at 
this point, the properties of the mixed valence state ap- 
pear to be very quite similar to those of the Kondo states. 

We also investigated the conductance properties of the 
SU(3) arrangement, proposed in Rcf. 9. At high temper- 
atures we observe in the side-conductance three charging 
peaks, corresponding to the three charging steps. As 
the temperature is lowered, the Coulomb blockade val- 
leys are gradually filled up, and two conductance shoul- 
ders emerge in the Kondo regime with a conductance 
G^(3/4)(e 2 /h). 

The methods and the computations presented here rep- 
resent a first and important step to perform DM-NRG 
and DMFT calculations for more elaborate SU(iV) or 
Sp(A^) lattice models. However, further work is neces- 
sary to optimize these DMFT calculations. While we 
definitely gain enormous storage space by using SU(iV) 
symmetries, the evaluation of the reduced matrix ele- 



12 



mcnts and multiple sums over internal representation la- 
bels are currently not carried out with maximal efficiency. 
Since for large JV's the dimensions of irreducible SU(A r ) 
representations grow very fast, these summations quickly 
become the bottleneck for DM-NRG calculations, and 
further work is needed to increase the efficiency of the 
corresponding subroutines. 
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N x N matrices 



Appendix A: Some details on SU(iV) representations 



In this appendix we give a brief overview of the rep- 
resentation theory of SU(iV), following the approach of 
Ref. 16, where a more detailed discussion can be found. 



While SU(2) has three generators, j z , j + , and j-, 
SU(iV) has iV 2 — 1 generators. We shall deal explicitly 
with only 3 (AT — 1) of them, denoted by Ji" and J± , 
where I = 1, . . . , N — 1. By definition, they satisfy the 
commutation relations 



i+i- 



2-> 

a+i- 



t+i- 



/0 
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/o 
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/o 



V 



i i+i 



\ o 
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I l+l 





1 



(A2a) 



0/ 
\ 



(A2b) 



0/ 



(A2c) 



0/ 



which generate the defining representation of STJ(iV). 
The N 2 — 3N + 2 remaining generators are obtained as 
commutators between generators involving different val- 
ues of I. Below we will use J a (and J a ) as collective 
notation for any of the generators (and corresponding 
matrices in the defining representation). 

Let |r, 7) denote the basis states of a general SU(iV) 
irrep, where T labels the irrep and 7 its individual basis 
states. The action of any generator on a basis state can 
be written as 



i a |r,7> = £ r 7',7liW), 



(A3) 



= ±j2\ 



= 2J«\ 



(Al) 



These have the same form as those of the corresponding 



where r°, = (r, 7'! J a \Y, 7) are the matrix elements 
of J a within the irrep T. For example, the matrix ele- 
ments of l ± in the defining representation are given by 
Eqs. (A2). 

It is convenient to choose the states |r, 7) to be si- 
multaneous eigenstates of ji^ for I = 1,...,N — 1, 
with eigenvalues Ap 7 , say. The sequence W z (r,j) = 

(X\ rt , ■ ■ ■ , A^/7_i) is called its weight. As for SU(2), all 
components of the weight takes half- integer values. 
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TABLE IV. Examples of GT-patterns and corresponding 
Young tableaux; the top row of each pattern determines the 
shape of its respective tableau. These examples have been 
constructed in such a way that each tableau/pattern con- 
tains the examples to its left as a sub-tableau/-pattern. Note 
that we usually drop columns of SU(JV) Young tableaux with 
length TV (e.g., the leftmost column in the rightmost tableau); 
this corresponds to subtracting the entry rriN,N from each en- 
try of a GT-pattern. The weight of a state can be directly 
constructed from a GT-pattern: Let <t; = Y2k=i i mfe J 
note the row sums, then Xf 1 — (ci+i — cr;)/2. 



FIG. 9. Weight diagram of the SU(3) irrep [p (non- 
orthogonal axes chosen to emphasize the symmetric structure 
of the irrep). Each dot represents a weight; we also indicate 
the Young tableaux of the corresponding states. The circled 
dot indicates a weight with inner multiplicity 2. The blue 
solid and dashed arrows represent the action of J_ and J_ , 
respectively. 



A convenient way to visualize all states of the same 
SU(TV) irrep is then given by weight diagrams (see Fig. 9 
for an example), which arc constructed by marking the 
point with coordinates W Z (T, 7) for each state |T, 7) of 
an irrep T. The operators J±* then map states onto their 



r,-> 



±1) 



neighbors in a weight diagram (A 

Each irrep T has a so-called highest-weight state 
(unique up to a phase), denoted by |r, 7 = T) for con- 
venience. It is annihilated by all Jl\ 



j2 ] 



|r, 7 = r) = o (z = i,...,tv-i). 



(A4) 



Its weight actually determines the properties of the whole 
irrep T, and is thus suitable to provide a labeling scheme 
for r. 

In contrast to SU(2), several states of an irrep T can 
have the same weight. The number of states with the 
same weight is called the inner multiplicity of this weight. 
Consequently, weights are not suitable as the label 7. 
Instead, we use one of two equivalent labeling schemes, 
Young tableaux or Gelfand-Tsetlin patterns. 

An SU(TV) Young tableau is a single, contiguous cluster 
of left-aligned boxes with at most TV rows, such that each 
row is not longer than the one above. Each box of a 
tableau carries a number between 1 and TV, inclusive, 
such that numbers do not decrease from left to right, 
and numbers increase strictly from top to bottom. 

A GT-pattern M is a triangular matrix of integers 



Mk i (1 < k < I < TV), commonly written as 



mijv-i • 



M 



mN-i,N-i 



\ 



(A5) 



/ 



which are subject to the so-called betweenness condition, 
mk,i > rrik,i-i > rrik+1,1 (1 <k <l < TV). (A6) 

Tab. IV gives examples of equivalent Young tableaux and 
GT-pattcrns. 

Young tableaux and GT-pattcrns arc composite labels, 
M = (T,7) in short, in the sense that they can play the 
role of both the irrep label Y and the state label 7. The 
shape of a Young tableau (i.e., without the labeling of 
boxes) determines an irrep F; it corresponds to the top 
row of a GT-pattern, whose A:-th entry rnk.N specifics the 
number of boxes in the fc-th row of the Young tableau. 
The dimension of an irrep Y is equal to the number of 
valid GT-patterns with a given top row. There exists a 
convenient formula for this number: 



dim(r) = 1 



l<k<k'<N 



mk.N - mk',N 



(A7) 



The GT labeling scheme has the advantage that for any 
of the generators J a G {J^i ±}, the corresponding matrix 
elements Y", [Eq. (A3)] within the irrep Y are known 
explicitly, given by a complicated formula worked out by 
Gelfand and Tsetlin. 32 

A further ingredient to SU(TV) representation theory is 
the decomposition of a tensor product Yi (g> T2 of two ir- 
reps into a direct sum of irreps (see Eq. (10)). For SU(TV), 
this trick is accomplished by the Littlewood- Richardson 
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rule, which is beyond the scope of this introduction, 
produces equations such as 



It 



EP®EP = D©EE|©2EP 



in 



(A8) 



The number of times a particular irrep occurs on the 
right-hand side is called its outer multiplicity; for SU(7V), 
it is > 1 in general. The particular basis transforma- 
tion effecting this decomposition is described by Clebsch- 
Gordan coefficients [see Eq. (11)]; Ref. 16 presents a nu- 
merical algorithm for computing them for any N . 

Now consider a quantum chain model involving N cre- 
ation and annihilation operators per site, /I™^ and 



for site n, with a — 1, 
the set of operators 



. , N. For a given site n, consider 



a,P=l 



(A9) 



with the matrices J a taken as the defining representa- 
tion of SU(JV) [Eq. (A2)]. The J a ^ satisfy the same 
commutation relations as the SU(iV) generators J a and 
hence generate an operator representation of SU(JV) on 
the Fock space of site n. The action of these operators 
partitions this Fock space into a direct sum of irreps. For 
example, Tab. V specifies these irreps explicitly for the 
case of SU(3). 

The Fock space of a the full chain is the direct prod- 
uct of the Fock spaces of each site. Correspondingly, 
J a = J a '^ generates an operator representation 

of SU(iV) on the, which can be decomposed into a di- 
rect sum of irreps by iterative use of the Littlewood- 
Richardson rule and Clebsch-Gordan coefficients. 

The Hamiltonian H for the full chain has SU(iV) sym- 
metry if it commutes with all generators J a . When the 
Hamiltonian is expressed in the Fock space basis just 
mentioned, consisting of a direct sum of SU(jV) irreps, 
it is block-diagonal, with each block containing matrix 
elements only between states transforming according to 
a given SU(iV) irrep. 

Diagonalizing such blocks, or more generally, calculat- 
ing matrix elements of operators, is expedited by using 
the Wigner-Eckart theorem. To this end, one needs to 
group operators in Fock space into operator multiplcts 
(sometimes called irreducible tensor operators). An op- 
erator multiplct transforming according to the irrep T is 
a set of operators Or, 7 that satisfy the relations 



J Q ,Or. 



7' 



(A10) 



For example, the set of all generators J a constructed in 
Eq. (A9) spans an operator multiplet, acting only on site 
n and transforming according to the adjunct representa- 
tion of SU(iV), which has dimension N 2 - 1. (For SU(3), 
this is the irrep [jp listed in Tab. I). 



Multiplets 


States 


Q F 


Energy 


Degeneracy 


1 


|1> = |o> 


-3/2 • 





1 


2 


|2)=4|0> ■ 
|3)=4|0> > 
|4)=4|0) . 


-1/2 □ 


£d 


3 


3 


|5)=44|0) " 

|6) = 44lo> > 

1 1) — a 2 «3l u / ) 


I — I 

1/2 □ 


2e d + U 


3 


4 


|8> = 444lo) 


1/2 • 


3e d + 3U 


1 



TABLE V. Organization of the impurity states into multi- 
plets. The 2 3 = 8 states are organized into 4 multiplets, each 
characterized by a set of quantum numbers (Q,F). States 
within a multiplet are degenerate in energy. 



A convenient way to explicitly construct an operator 
multiplct associated with a given site n is to find its 
highest-weight operator by guessing or repeated appli- 
cation of raising operators, and then produce the other 
operators in the multiplet by applying lowering opera- 
tors, using Eq. (A10). Tab. I gives some examples of 
SU(3) operator multiplets constructed in this manner. 

The tensor product of two operator multiplets, each 
acting on individually on a separate site, can be de- 
composed into a direct sum of two-site operator mul- 
tiplets, again using the Richardson-Littlewood rule and 
Clebsch-Gordan coefficients. However, one often finds 
that extending a single-site operator multiplct to two 
sites by taking its tensor product with the identity op- 
erator (transforming according to the trivial SU(iV) rep- 
resentation) is enough, so a complicated decomposition 
can be avoided in most cases. For more than two sites, 
this procedure is applied iteratively. 



Appendix B: Multiplets within the Anderson model 
with Uq(1) x SUf{3) symmetry 

In this appendix we illustrate in more detail how the 
general concepts presented in Appendix A can be applied 
to our case, when Uq(1) x SUf(3) symmetry is used. We 
first construct the lowering/raising operators in the Fock 
space by using (A9). Explicitly, for SU(3) we shall need 
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explicitly six generators, 

J_l_ ' = d\d,2 1 j\_ ^ = G?2<^3j 

Ji 1] = \ (4di - d\d 2 ) , jC 2 > = 1 (4d 2 - djda) . 

The initial impurity states can be constructed and or- 
ganized into 4 SU(3) multiplets relatively easily. They 
are presented in Table V. The highest weight states 
|r,7 = r) can be found by requiring that both Jv' and 

"(2) 

J_j_ annihilate them. Further states within the multi- 

plct arc then obtained by acting with and and 
comparing to Eq. A3. Notice that the Clcbsch-Gordan 
coefficients do depend on the particular choice of basis, 
and therefore, to make use of the Wigner-Eckart theorem, 
every multiplct must be constructed to conform with the 
same choice of basis as the Clebsch-Gordan coefficients. 
This choice of basis is implicitly contained in the ma- 
trix elements (I± ) 7 y. In the open access Flexible NRG 
code 13 we used the conventions of Ref. 16, for which these 
matrix elements are explicitly given in Ref. 32. 

To use the Wigner-Eckart theorem, we need to orga- 
nize operators into operator multiplets. A few examples 
of these were presented in Tab. I. Let us discuss here the 
specific case of the spin operator, forming a 8-dimcnsional 
multiplct, which transforms according to the [p repre- 
sentation. A similar procedure applies to the hopping or 
the superconducting operators. 



We can quickly guess that the highest weight opera- 
tor is d\d^ since it commutes with both J_l and 
(see Eq. (A10)). This operator thus corresponds to the 
[|P state. Now, by applying the lowering operators (i.e., 
forming commutators with them) as indicated in Fig. 9, 
we can generate all the other operators that form the 
multiplct. 16 In this procedure one needs again the explicit 
form of the corresponding matrix elements, (I±) 77 / to 
fix the proper phase/sign of the basis states and their 
normalization. 

We remark that particular care must be taken in cases 
where a weight has an inner multiplicity, as for the mem- 
bers [|P and [§p of the spin operator multiplet. There one 
needs to take the correct linear combination of the cor- 
responding operators to transform according to (I± ) 77 / 
(see Eq. (A10)). 



Appendix C: Recursion formulas 

In this section we shall detail how the recursive rela- 
tions (Eqs. (16)-(19)) were derived. The general proce- 
dure is based on the Wigner-Eckart theorem and using 
sum-rules satisfied by the Clebsch-Gordan coefficients. 
We first derive the recursion relation for the irreducible 
matrix element of the hopping Hamiltonian. The general 
expression of the hopping matrix element at iteration n 
reads 



{u,r u ,^/ u \r n -i tn \v,T v ,j v )^ = y 



(Cl) 



Here the multiplets, \u,T u ,7 u y™' and |u,r„,7„)^ were 

constructed in terms of the block multiplet \i, r i; ji)^ n ^ 
and local multiplet rj^ c , 7^°°) at iteration n — 1 using 



I 

Eq. (11). Next we exploit the locality of the C™\ , i.e. 
that the operator C^™ r acts only on local states 
while (C^jT 1 ^ y acts on the block states i. The matrix 
element of Eq. (Cl) then becomes 



,r„,7 U | r„_i, n |u,r„,7„)^| = E 



sgn(C a ,fx) E E (rJr.T^ri.T* I r -7n)l u x (C2) 



(u, r>° c , 7 loc | eft | M , , 7 >° c ) + h.c. 



Jn-l] 



Here the sign function sgn(C a , /i) = ±1 arises as we com- 
mute the local state fj, over the operator Cl™ -1 !. If the 
hopping operator C is fermionic and the local state con- 



tains an odd number of fermions, then the sign is nega- 
tive, otherwise it is positive. Now we can use the Wigner- 
Eckart theorem, Eq. (13), and express the matrix ele- 



16 



merits of the creation/annihilation operators in terms of their reduced matrix elements. By doing that, we imme- 
diately recover the result (16) with 



D(a,a,/3;u,v) = sgn(a, M )E E E (1^,7^,7* I r„,7„)^ (r 1 -, 7 '° c ; T,, 7 , | T u , lu )^ x 

7i.7^ oc 7j,7l° c 

(r B) 7a ; r is7j I r,, 7 )I (r„, 7 a ; rjr, 7 jT I rj? c , 7'°% . (C3) 

A slightly different analysis can be done for block operators. First we notice that we can "invert" the Wigncr-Eckart 
theorem, Eq. (13), using the completeness of the Clebsch-Gordan coefficients, and express the reduced matrix elements 
instead as: 

(u i| A || v)f ] = E I a taha k r v,iv)al ( r A, 7a; r„, 7 „ I v u , J u )p . (C4) 

1A, ~1v 



As a next step, we need to express the matrix element 
7u| ^r. 4 ,7A l u ir c ,7t))^ in (C4) in terms of the 
irreducible matrix elements of the operator at iteration 
71—1. Here we follow the same strategy as in the case 



of the hopping operator: We expand first the states us- 
ing (D2), while keeping in mind that the operator acts 
only in the block states sector, and then use the Wigner- 
Eckart theorem (13) for the operator's matrix elements. 
By doing that we find: 



a„ (u, r u , 7„| Ar A , yA \v, T v , 7„)W = E E ^."^c^joc^oc^ioc (r}° c , 7^ oc ; T l7 lt \ T u , 7„)*^ 

7i , 7j 7 loc )7 loc 

(^,7^5^,7, I r„, 7 «) Q „ >< 

sgn(A/i)E( r ^'^ r ^^' I II A II it~ l] ■ (C5) 

a 

Plugging Eq. (C5) in Eq. (C4) we recover the result for the block operator stated in (18) with the coefficient F defined 
as 

f (a, /?;«,«) = 8gn(A,/i) E E E ( r lT> 7^ c ; r 4) 7, 1 r„,7«)^ (rjr, 7jT; r„ 7j I r„,7„) a . x 

7 loc 7i , 7:j 7A , 7 j 

(rA,7A;rj,7j | r i)7i )* (i\4,7A;r„, 7 „ I r^.^. (C6) 

A similar analysis can be done in the case of a "local operator" , giving the final expression for the coefficient K 
entering Eq. (19) 

K(a,(3;u,v) = J2 E E (^° c , 7^; T,, 7 I T u , 7^^ (r^ c , 7 ^ c ; T,, 7 | T w , 7 „)^ x 

(r Aj 7^1^,7^ I rjr,7j? c )* (1^,7^,7, 1 r„, 7 „V (C7) 



Appendix D: Connection to the Matrix Product 
States approach 

It is well-known that the states constructed within the 
NRG framework can be viewed as Matrix Product states 
(MPS) 33 . Moreover, it has been shown recently that non- 
Abelian symmetries can be incorporated into the con- 
struction of MPS 15,34 . In this appendix we briefly review 
how this can be done in the context of NRG. We start 



with the simple observation that the Hamiltonian of a 
chain of length N, Hn acts on the Hilbert space spanned 
by a basis constructed from local multiplets: 

{1^0,1^,7^0 , Mi, 1 Ml ,7 M1 , ■ • ■ MN ,7 MJV /; • 

(Dl) 

The dimension of this basis set is dod N , with do the di- 
mension of T~Lq, and d is the dimension of the Hilbert 
space at sites along the Wilson chain (see Tab. V). For 
the SU(3) Anderson model we have four multiplets on 
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each site, and do — d = 8. 

Let us now assume that we have constructed somehow 
some block states i, which span the relevant part of the 
Hilbcrt space of a chain with n — 1 sites. The number 
of these states, D is, of course, much less than the total 
number of states within this block, which would be of the 
order of ~ d n . We can then use these states to express 
the appropriate (relevant) states for a chain of n sites as 



E 

M,rri, T r ] 



p 



r\7^,r^7^;o:- 



l/*,r)r,7jr>®|i.ri,7i> [n - 1] • (D2) 

Here P is some projector that generates the relevant 
block of multiplcts from the block multiplcts {\i)} 
by adding some local states The tilde in the sum 

indicates that only a number D of multiplets are kept at 
each iteration. In the presence of symmetries, the matri- 



ces P can be factorized into products of reduced matrix 
elements and Clcbsch-Gordan coefficients, as: 



rploc 



fiTi 



with 



lill 



(D3) 
(D4) 



Here the key observation is that the reduced matrix 

rp loc 1 r 1 

elements (^4 r .p_ ) can also be thought of as the matrix 

elements of some irreducible operators, labeled by T l ° c . 

Repeating the iteration procedure multiplcts describ- 
ing the full chain of size N can be constructed this way as 
a matrix product state (summation of repeated indices is 
implicitly assumed) 



IV) [JV] = lv,r^ I 7^> 



[N] 



rploc loc, 
J^ L Mo ' 'CD J 



[mo] 



A 



rploc. 



ll-'l 



c 



rploc. 



7>n7i; 



.4 



ra 



ra 



Tii7« : 



x ... X ( A [r ^ 



r, r, 

*iV-l % 



rploc , 
Q l MN J 

^'JV — lT*^ ' ai i/; 



f7 1 
I 'MiV J 



7i W „ 1 7i^, 



I,, t^Ioc _,loc. ,, t^Ioc _,loc. . ,, jnloc „loc \ 

lMo,l Mo ,7 Mo jMi, 1 ^ ,7^ >■■•■> (*n, L m>1hn) 



Note that the index structure that arises here implies 
matrix multiplication not only for the A-matrices of re- 
duced matrix elements, but also for the C-matrices of 
Clebsch-Gordan coefficients. This MPS formulation can 
be used either to implement standard Wilsonian trunca- 
tion (as done in Ref. 15), or, alternatively, to proceed 
variationally, as done in the density-matrix renormaliza- 
tion group (DMRG), 35 whose use for Wilson chains was 
explored in 33 ' 36 . In the latter case, one views the A ma- 
trices as a set of variational parameters that need to be 
optimized according to some criteria. To find the optimal 
approximation for the ground state, e.g., we look for the 
corresponding MPS which minimizes the total energy, 



E 



This variational problem can be converted into a gener- 
alized eigenvalue problem and solved using an iterative 
sweeping-like procedure (see [35] for more technical de- 
tails). Once the MPS state is found, it is possible to 
construct the eigenspectrum of the Hamiltonian ipv for 
a fixed N. The flow diagram, such as the one presented 
in Fig. 2, can be obtained by tracing the spectrum of the 
Hamiltonian with increasing N; if the discretization pa- 
rameter A is large enough (A > 2), the numerical results 
thus obtained are essentially equivalent to those using 
standard Wilsonian truncation 36 . For more details on 
how the operators can be treated at the MPS level we 
direct the reader to the more detailed reviews, Ref. 15 
and 35. 
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